Circuit QED scheme for realization of the Lipkin-Meshkov-Glick 
model 
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Abstract. - We propose a scheme in which the Lipkin-Meshkov-Glick model is realized within a 
circuit QED system. An array of A'' superconducting qubits interacts with a driven cavity mode. 
In the dispersive regime, the cavity mode is adiabatically eliminated generating an effective model 
for the qubits alone. The characteristic long-range order of the Lipkin-Meshkov-Glick model is 
here mediated by the cavity field. For a closed qubit system, the inherent second order phase 
transition of the qubits is reflected in the intensity of the output cavity field. In the broken 
symmetry phase, the many-body ground state is highly entangled. Relaxation of the qubits is 
analyzed within a mean-field treatment. The second order phase transition is lost, while new 
bistable regimes occur. 



Introduction. — During the last two decades, sys- 
tems involving single atoms or superconducting qubits in- 
teracting with a quantized mode of a high-Q resonator 
have supported some of the finest experiments analyz- 
ing fundamental aspect of quantum mechanics. Among 
these cavity /circuit quantum electrodynamics (QED) ex- 
periments, it is especially worth mentioning; entanglement 
generation fTlE] , quantum-classical transition 131 , and veri- 
fication of the graininess of the quantized electromagnetic 
field [3]. Very recently, a new branch of cavity /circuit 
QED including several atoms/qubits has emerged. Co- 
herent coupling of Bose-Einstein condensate [5] as well 
as coupling of few qubits [SIISIIZ] to a single cavity mode 
have been experimentally demonstrated. These experi- 
ments paves the way for controlled studies of novel coop- 
erative phenomena of many-body systems. 

Microscopic many-body systems are at the heart of a 
wide range of physical branches, e.g. condensed matter, 
quantum information processing, Bose-Einstein condensa- 
tions, and ultracold atoms. With experimental progress, 
much attention has been especially paid to phase tran- 
sitions and more recently to ground state entanglement 
properties. At zero temperature, a phase transition in 
these systems is driven by quantum fluctuations and is 
thereby termed quantum phase transition [5]. In the realm 
of quantum information processing, determining the quan- 
tum properties of the many-body ground state is of cen- 



tral significance and especially the amount of entangle- 
ment distributed among the entities 9,. Spin models, 
describing interaction between spin 1/2 particles are of 
special importance as each spin can be ascribed a qubit 
state. Most common are lattice models, e.g. Ising and 
XX models ^8^ , characterized by short or semi-short range 
interaction. Lately, also long range models like the Lipkin- 
Meshkov-Glick (LMG) have gained extra recognition since 
they might be realizable in the near future in various sys- 
tems. 

The LMG model has its origin in nuclear physics [TO] . 
In the recent past, several proposals for realization of the 
LMG model in different contexts have been put forward, 
for example within trapped ions JJJ, Bose-Einstein con- 
densates in double- wells [T5] or in cavities [T3], or of non- 
interacting atoms / qubits within a cavity |14lfT5] . The gen- 
eral properties of the LMG model have been considered 
by numerous authors, focusing on entanglement proper- 
ties [12] and finite size effects |17) . 

In this paper we show how the symmetric LMG model 
may be achieved by extending the recent experiment of 
ref . [5] ; an array of superconducting qubits are identically 
coupled to a cavity mode. The resonator is externally 
pumped, and in the dispersive regime the cavity field is 
adiabatically eliminated yielding an effective model for the 
qubits. The quantized cavity mode serves as a quantum 
bus |2j, which induces the long range interaction between 
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the qubits. For zero couphng of the qubits to its environ- 
ment, we consider the full many-body problem and cal- 
culate the ground state entanglement both for the LMG 
model and an extended model also realizable in the same 
setup. Decoherence of the qubits is analyzed within a 
mean-field approach. We demonstrate that bistable dy- 
namics emerges in this situation, while the second order 
quantum phase transition of the LMG model is lost. All 
these ground-state properties are accessible via detection 
of the output cavity field leaking out of the end mirrors. 

Model system. — Inspired by the recent circuit QED 
experiments and the quest for engineering multi- 

qubit systems, we consider an array of N superconducting 
qubits coupled to a single mode of an optical resonator. 
The individual effective qubit-field coupling strength gi 
[i = 1, 2, ...,N) are taken to be all equal, i.e. gi = g, 
V i. This is in agreement with ref. [6] where the fluctu- 
ations of the (7i's were only a few percent. The quality 
factor Q of the resonator is assumed large enough to sup- 
port well resolved mode frequencies, but still not infinite 
such that the photon decay rate k 7^ 0. The cavity is 
externally driven via one side mirror by a pump with fre- 
quency LOp and amplitude 77. Considering an interaction 
frame rotating with the pump frequency and introducing 

u!p and the atom- 



the cavity-pump detuning Ac — uJc 



pump detuning Aa 



the Hamiltonian is written 



H 
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Ho + Hc 
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cav 



Here, a and are annihilation and creation boson oper- 
ators of the cavity mode respectively, = J2iLi 
is the collective spin operator of component a {— x, y, z) 
with ai*^ = \e)ii{g\ + \g)ii{e\ and ai*^ = \e)ii{e\ - \g)ii{g\ 
the regular Pauli matrices acting on the ground \g) and 
excited |e) state of qubit i, and finally Hcav and Hqu repre- 
sent coupling to surrounding reservoirs leading to respec- 
tively photon leakage and qubit decoherence. 

In the absence of losses and pumping, is the Dicke 
Hamiltonian introduced to describe superradiance of N 
two-level atoms interacting with a quantized cavity mode. 
If g < ^uja^^c, the ground state of the system is said 
to be in the normal phase charactericed by vacuum cav- 
ity field and all qubits in their ground states, while if 
g > y'ujaUJc the ground state is an entangled superradi- 
ance phase where the field is in a Schrodinger cat and the 
qubits partly excited [18]. The quantum phase transition 
between the two phases is of second order nature without 
the rotating wave approximation. The qubit-field coupling 
is inversely proportional to the square root of the effective 
mode volume V, implying that if the density po — N/V 
is fixed we may write g = g^/^fN. The scaling of the 
coupling gives a well defined thermodynamic limit. In the 



presence of external pumping, one must assume 77 — 770 VjV 
in order to have a proper thermodynamic limit. Note that 
this is intuitive, if the cavity volume is enlarged the pump- 
ing should be increased accordingly. 

If qubit relaxation 7 and qubit dissipation 7' can be set 
to zero within the time-scales of interest, the system state 
p obeys the master equation 
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Heisenberg equations of motion become 

a = -(k -I- iAc)d + rjoVN - i-^^Sx + V2Kam{t), 

V N 

Sx = -AaSy + (d + an Sz, 



N 

(3) 

Sy = AaSx + i—^ (a - d^) 

Sz = (a - a^) Sy - (a + a^) Sx- 

Here, dot represent time-derivative and di„(i) is the in- 
put Langevin noise obeying (ai„(i)) = (d-„(t)di„(i')) = 
and (dm(t)d|'„(i')) = S{t - t') [19]. In what follows, the 
Langevin force din {t) will be omitted. In order to derive an 
effective spin Hamiltonian we assume a separation of time- 
scales, the cavity field is taken to evolve on a much shorter 
time scale than the qubits. This implies \k + iAc\ 3> go, 
and in this regime we may consider the steady state of the 
cavity field 
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Inserting this expression, together with its corresponding 
hermite conjugate, into the equations of motion for the 
collective spin operators the resulting equations defines 
an effective Hamiltonian for the qubits. To order iV" one 
obtains 



^eff = hAaSz - hhSx 



(£,2 



St 



where 



h = - 



2gor]oK 



A = 



(5) 



(6) 



For driving at resonance with the atomic transi- 
tion, Aa = 0, the above Hamiltonian is unitary 
equivalent with the isotropic LMG model via U = 



exp 



The long-range cou- 



f (5, + Sy + 5, j /V3 
pling, inherent in the squares of the collective spin op- 
erators, is a direct outcome of the cooperative qubit-field 
interaction. The second term, proportional to Sx , charac- 
terizes driving of the qubits, and derives from the driving 
of the cavity field. Such a term is expected since driving 
of the qubits or field is unitary equivalent in the Dicke 
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model. We note that a special case of the LMG model has 
been discussed in the context of circuit QED earlier, with 
the focus on few qubits and a closed system [15] and not 
on phase transitions nor bistability. 

Many-body properties of the effective qubit sys- 
tem. — The isotropic LMG model, i.e. Aa = in 
eq. ([5]), is analytically solvable. Since the total spin 
commutes with Heff and so does Sx, the eigenstates of 
Heff are \S,M)x, where S^\S,M)x = 5(5 + l)|5,Af)^ 



and Sx\S,M) = M|5, Af)^. 
cle 



The eigen energies per parti- 
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where the ground state is to be found in the maximum 
spin sector, 5 = N/2, and 
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\h\ < A 
\h\ > A. 



(8) 



Here, I{x) is the integer part of x rounded off to nearest 
integer. Without loss of generality, we have restricted the 
analysis to the ferromagnetic case A > 0. (Sx) — Mq being 
an order parameter, the ground state possess a second 
order phase transition at he = ±A. For \h\ > A the ground 
state is ferromagnetic, while in the broken symmetry phase 
\h\ < A the system exhibits an multi-partite entangled 
ground state. 

For a symmetric multi-partite qubit state, the reduced 
density matrix p2q for two of the qubits is easily obtain- 
able |20j . The concurrence is a proper measure of en- 
tanglement for a general, mixed or pure, two qubit state. 
Given p2q, we introduce the operator g2q = P2q{o'iy ® 
^'^y) P2qi^iv ^2y) where aiy is the Pauli y-matrix acting 
on qubit i and asterix denotes complex conjugation. The 
concurrence is defined as 



C = max {0, fii - p2- - M4} , 



(9) 



with /ij the square root of the j'th eigenvalue of g2q or- 
dered such that /ii > /i2 > Ms ^ M4- For ^ disentangled 
state C = 0, while C = 1 characterices a maximally entan- 
gled one. In the case of the isotropic LMG model one has 
the rescaled concurrence, Cr — {N + 1)C, given by [20] 



Cr = 



1 

2N 



^{N^ - 4Af 2)[(7v - 2)2 - 4M2] I . 



(10) 



For our numerical results we try to employ physically 
relevant parameter values. Experimentally, the photon 
decay rate k is of the order 1 to 50 MHz OEl]. The 
qubit-field coupling ~ 600 MHz in the strong coupling 
regime experiments [SI HI]- In the following we will ex- 
press frequencies in units of k, and take go — 100k. Fur- 
thermore, to assure validity of the adiabatic elimination 
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Fig. 1: (Colour online) The order parameter ~ {Sx)/N (a) 
and scaled concurrence Cr (b). For the red line Aa = k/5, 
while for the black one Aa — 0. It follows that the phase tran- 
sition at \h\ — \c = 5k is lost once Aa 7^ 0. The parameters 
are go = 100k and Ac — 2000k. 



of the cavity field, we use Ac = 2000k, i.e. 20 times 
larger than the qubit-field coupling. In fig. [U the or- 
der parameter Sx = {Sx)/N and scaled concurrence Cr 
are shown as function of h/n. Varying of h is preferably 
achieved by tuning the pump amplitude 770 since /i cx 779. 
In the two plots, the black lines depict results for the 
LMG, i.e. Aa = 0, while the red lines give the same for 
Aa = k/5. The number N of qubits is here 200, but the 
results are relatively insensitive to it. That is, even for 
few qubits is the scheme efficient for entanglement gener- 
ation. We note that the second order phase transition at 
1^1 = Ac = 5k vanishes for non-zero A^. For the LMG 
model, the rescaled concurrence is approximately 1 in the 
broken symmetry phase indicating maximal distribution 
of the entanglement between the qubits. In the vicinity 
of the critical points, the concurrence increases abruptly, 
hence giving a clear signal of the phase transition. The 
amount of quantum correlation for \h\ < Ac is lower when 
Aa 7^ 0, but contrary to the LMG model, the concur- 
rence survives also for regimes with \h\ > Ac. However, at 
h sa 6.5k the entanglement suddenly vanishes. This rapid 
disappearance is only present in the concurrence and not 
in the order parameter Sx , and is therefore not an outcome 
of a phase transition. Indeed, sudden death of entangle- 
ment has been studied in great detail lately [22^ . 

Apart from constituting an effective qubit-qubit inter- 
action, the quantum properties of the cavity field allow for 
non-destructive detection of the multi-partite qubit state. 
The amount of leakage of photons out of the end mirror 
is directly proportional to the photon number in the res- 
onator ^23j. Since, 
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measuring n^s gives a handle of the qubit state. In partic- 
ular, in the LMG case we have (5^) — Mq for the ground 
state with Mq given in eq. (|8]). 

Influence of qubit decoherence. In the case of 
two-photon Raman coupling via a highly detuned third 
level, the two qubit states \g) and |e) can be made insensi- 
tive to decay. However, for a direct transition, the excited 
|e) state will decay with some rate 7. Such irreversible 
processes will affect the system properties. The collective 
decoherence and dissipation are modeled with a master 
equation for the qubits, 

+ {2S+PqS- - S-S+Pq - PqS-S+^ 




where S± — i Sy and 7' characterizes the decoherence 
rate. Hence, the 7-term describes decay of the excited 
qubit level and the 7'-term loss of coherence due to inter- 
action with its environment. We will assume the thermo- 
dynamic limit, where we can neglect quantum fluctuations 
and thereby factorize (AB) = (A) (B) for any operators A 
and B. Hence, we consider a mean- field approach. De- 
noting Sa = {Sa)/N, a — X, y, z, we obtain the following 
semiclassical equations of motion 

Sy ^ AaS^ + hsz - Xs^s^ - jSzSy, (13) 

Sz = -hSy + XSySx +l{sl + 4) 

in the thermodynamic limit. In ()13p . only terms of order 
are kept, and as a consequence contribution from the 
dissipation vanishes. Throughout this section we will set 
Aq — 0. In addition to (IT51) . we have the angular mo- 
mentum conservation constrain + Sy + = 1. The 
system of equations (|13|) posses in general several fixed 
points {sx,Sy,Sz). Analyzing the fixed points, i.e. study 
the steady state solutions {sx — Sy — = 0) of ([T3l) . 
one finds that if 7 = the many-body results of the pre- 
vious section are regained. If, on the other hand, 7^0 
it directly follows that either Sx or Sz must be identically 
zero. 

For Sx — 0, the steady state solution is easily found as 




valid for \h\ < 7. In the other case, — 0, the steady 
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Fig. 2: The steady state solutions of the truncated Heisen- 
berg equation (|15|) . The solid lines are stable solutions, while 
dotted lines represent unstable solutions. The parameters are 
go = 100ft:, Ac = 2000k, and 7 = k/5 (a) or 7 = ft:/50 (b). 

state solutions obey 




This set of solutions shows an imperfect transcritical bi- 
furcation pattern, see [2l]. The last constant term (7/A)^ 
cause the bifurcation to become imperfect. As an out- 
come, the phase transition encountered in the lossless case 
(7 = 0) at = Ac is lost and the transition is instead 
smooth. In addition, multiple solutions may exist of which 
some are stable and others unstable. The stability is found 
via a linear stability analysis around the fixed points |24] . 
i.e. perturb around the steady state solutions and check 
whether the perturbations increase or decrease with time. 
In fig.[2]we display the steady state solution Sx correspond- 
ing to the solutions (fT5|) for 7 = k/5 in (a) and 7 = k/50 
in (b) . The value of 7 in (a) gives a realistic decay rate of 
an excited state in a qubit [21]. In the figure, solid lines 
represent stable solutions, while the unstable solutions are 
marked by dotted lines. 

It is seen in fig. [5] that a maximum of three stable so- 
lutions may coexist given a set of parameters. The one 
solution picked by the system is determined by details of 
the preparation and dynamical process taken by the ex- 
periment at hand. For example, by dynamically varying a 
system parameter the evolution must not imply that the 
system follows the ground state. However, if the charac- 
teristic time scale of the evolution is long compared to the 
inherent time scale for quantum fluctuations it is likely 
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that transitions from stable solutions with higher energy 
will decay into stable solutions with a lower energy. Of 
the three stable solutions of fig. [2l the solution with small- 
est \sx\ is presumably having the lowest energy. This as- 
sumption derives from estimating the energy utilizing the 
energy density £ = —hhsx — h\ (s^ -I- s^). Thus, the cou- 
pling of the qubit to the environment enters in the actual 
value of the mean- field values (s^, Sy,Sz) and not directly 
in the energy density. When the solid lines terminate at 
\h\ w K of fig. [2] (a), the corresponding solutions become 
complex. It therefore follows that there may occur pa- 
rameter regions in which no stable solution exists. The 
solution ([T4)) . valid for < 7, is stable but still its exis- 
tence does not warrant stable solutions within all param- 
eter regimes. Thus, tuning h across k from \h\ > k, 
and assuming the system to follow adiabatically its in- 
stantaneous ground state, will cause sudden jumps in the 
physical quantities. This bistability will be reminiscent 
of a first order dynamical phase transition. It should be 
pointed out that these results are derived in the thermo- 
dynamic limit, hence it does not necessarily imply that 
the finite size system also lacks stable solutions within 
these parameters. As discussed in the previous section, 
the output cavity field being proportional to n^s contains 
information about the actual state of the qubits. This, of 
course, holds also in this mean-field treatment and thereby 
one should be able to detect the dynamical bistability in- 
dicated by fig. m The second (b) plot makes clear that 
in the limit of no qubit decay, 7 = 0, the results of the 
previous section is captured. We note that optical bista- 
bility has gained renewed interest lately in optomechanical 
systems pS] . 

Conclusion. — We have shown that under proper con- 
ditions, the LMG many-body spin model can be realized 
within circuit QED by considering an array of supercon- 
ducting qubits coupled to one quantized mode of a higli-Q 
resonator. The symmetric LMG model is obtained when 
the cavity pump is resonant with the qubit transition, 
and it was demonstrated how the cavity output field re- 
veals sufficient information for determining the phase of 
the many-body qubit ground state, i.e. the ferromagnetic 
or the broken symmetric one. The many-body entangle- 
ment properties of the broken symmetric phase was dis- 
cussed. Driving the resonator off resonant with respect 
to both the cavity and the qubits, the phase transition 
of the LMG model is lost and the ground state becomes 
in general entangled also in the otherwise ferromagnetic 
phase. 

Within a mean-field treatment we considered the sit- 
uation where qubit relaxation and dephasing are taken 
into account. An imperfect bifurcation structure emerges 
where the second order phase transition of the LMG model 
vanishes, but bistable properties appear. Driving the sys- 
tem through such bistable points would probably lead to 
sudden changes in the measured quantities. 

In our treatment, all qubits were assumed identical and 



also to couple identically to the cavity mode. In ref. [B], 
an experiment containing three qubits was manufactured 
and it was found that the fluctuations of the couplings 
were only a few percent. For larger number of qubits and 
greater fluctuations, new phenomena is likely to occur. Of 
particular interest is the realization of spin glasses. Spin 
glasses are normally outcomes of disorder within spin mod- 
els [26]. Thus, the considered model might be suitable for 
controlled generation of spin glass systems. An advantage 
of considering superconducting qubits rather than ultra- 
cold atoms is the absence of particle motion. It is known 
that the additional degrees of freedom induced by parti- 
cle motion can greatly affect the system dynamics and its 
phase diagram properties [27]. 

* * * 

The author acknowledges financial support from the 
Swedish research council- Vetenskapsradet/VR. 

REFERENCES 

[1] HAGLEY E., MAITRE X., NOGUES G., WUNDERLICH C, 

Brune M., Raimond J.-M., and Haroche S. , Phys. 
Rev. Lett, 79 (1997) 1; DiCarlo L., Chow J. M., Gam- 
betta J. M., Bishop L. S., Johnson B. R., Schuster 
D. I., Majer J., Blais a., Frunzio L., Girvin S. M., 
and Schoelkopf R. J., Nature, 460 (2009) 240. 

[2] Majer J., Chow J. M. , Gambetta J. M., Koch J., 
Johnson B. R., Schreier J. A., Frunzio L., Schuster 
D. I., HouCK A. A., Wallraff A., Blais A., Devoret 
M. H., Girvin S. M., and Schoelkopf R. J., Nature, 
449 (2007) 443. 

[3] Brune M., Hagley E., Dreyer J., Maitre X., Maali 

A., WUNDERLICH C, RAIMOND J.-M., AND HAROCHE S. 

, Phys. Rev. Lett, 77 (1996) 4887 

[4] Rempe G., Walther H., and Klein N. , Phys. Rev. Lett., 
58 (1987) 353; Brune M., Schmidt-Kaler F., Maali A., 
Dreyer J., Hagley ., Raimond J.-M., and Haroche S. 
, Phys. Rev. Lett, 76 (1996) 1800; Schuster D. I., Houck 
A. A., Schreier J. A., Wallraff A., Gambetta J. 
M., Blais A., Frunzio L., Majer J., Johnson B., De- 
voret M. H., Girvin S. M., and Schoelkopf R. J., 
Nature, 445 (2007) 515; Fink J. M., Goppl M., Baur 
M., BiANCHETTi R., Leek P. J., Blais A., and Wall- 
raff A., Nature, 454 (2008) 315. 

[5] COLOMBE Y., Steinmetz T., Dubois G., Linke F., 
Hunger D., and Reichel J. , Nature, 450 (2007) 272; 
Brennecke F., Ritter S., Donner T., and Esslinger 
T. , Science, 322 (2008) 235. 

[6] Fink J. M., Bianchetti R., Baur M., Goppl M., Stef- 
FEN L., FiLipp S., Leek P. J., Blais A., and Wallraf 
A., Phys. Rev. Lett, 103 (2009) 083601. 

[7] Altomare F., Park J. L, Cicak K., Sillenpaa M. A., 
Allman M. S., Li D., Sirois A., Strong J. A., Whit- 
taker J. D., AND SiMMONDS R. W., arXiv:1004.0026, (.) 

[8] Sachdev S., Quantum Phase Transitions (Cambridge Uni- 
versity Press, Cambridge) 1999. 

[9] Amico L., Fazio R., Osterloh A., and Vedral V., 
Rev. Mod. Phys., 80 (2008) 517. 



p-5 



Jonas Larson 



[10] LiPKiN H. J., Meshkov N., and Click A. J., Nud. 

Phys., 62 (1965) 188. 
[11] Unanyan R. G. and Fleischhauer M., Phys. Rev. Lett., 

90 (2003) 133601. 
[12] MiCHELi A., Jaksch, Cirac J. I., and Zoller p., Phys. 

Rev. A, 67 (2003) 013607. 
[13] Chen C, Liang J.-Q., and Jia S., Opt. Express, 17 

(2009) 19682. 

[14] Morrison S. and Parkins A. S., Phys. Rev. Lett., 100 
(2008) 040403; Phys. Rev. A, 77 (2008) 043810. 

[15] TSOMOKOS D. I., ASHHAB S., AND NoRi F., New J. 
Phys., 10 (2008) 113020. 

[16] ViDAL J., Palacios C, AND AsLANGUL C, Phys. Rev. 
A, 70 (2004) 062304; Latorre J. I., Orus R., Rico E., 
AND ViDAL J., Phys. Rev. A, 71 (2005) 064101; Orus R,., 
DUSUEL S., AND ViDAL J., Phys. Rev. Lett, 101 (2008) 
025701; Cui H. T., Phys. Rev. A, 77 (2008) 052105. 

[17] DusuEL S. AND ViDAL J., Phys. Rev. Lett, 93 (2004) 
237204; Phys. Rev. A, 71 (2005) 224420; KwOK H. M., 
NiNG W. Q., Cu S. J., AND Lin H. Q., Phys. Rev. E, 78 
(2008) 032103. 

[18] Wang Y. K. and Hioe F. T., Phys. Rev. A, 7 (1973) 
831; Hepp K. and Lieb E. H., Phys. Rev. A, 8 (1973) 
2517. 

[19] Cardiner C. W. and Zoller P., Quantum Noise 
(Springer, Heidelberg) 2004. 

[20] Wang X. and M0ller, Euro. Phys. J. D, 18 (2002) 385. 

[21] Schuster D. L, Houck A. A., Schreier J. A., Wall- 
raff A., Cambetta J. M. , Blais A., Frunzio L., 
Johnson B., Devoret M. H., Cirvin S. M., and 
Schoelkopf R. J., Nature, 445 (2007) 515. 

[22] Yu T. and Eberly J. H., Phys. Rev. Lett, 93 (2004) 
140404. 

[23] Boozer A. D., Boca A., Miller R., Northup T. E., 
and Kimble H. J., Phys. Rev. Lett, 97 (2006) 083602; 
FoRTiER K. M., Kim S. Y., Cibbons M. J., Ahmadi P., 
AND Chapman M. S., Phys. Rev. Lett, 98 (2007) 233601. 

[24] Strogatz S. H., Nonlinear Dynamics and Chaos (West- 
view Press) 1994. 

[25] Larson J., Damski B., Morigi C, and Lewenstein 
M., Phys. Rev. Lett, 100 (2008) 050401; Jayich A. M., 
Sankey J. C, ZwiCKL B. M., Yang C, Thompson J. 
D., Cirvin S. M., Clerk A. A., Marquardt F., Har- 
ris J. C. E., New J. Phys., 10 (2008) 095008; TiAN L., 
arXiv: 1003.2006, (.) 

[26] Jarrett T. C, Lee C. F., and Johnson N. F., Phys. 
Rev. B, 74 (2006) 121301. 

[27] Larson J. and Lewenstein M., New J. Phys., 11 (2009) 
063027. 



p-6 



